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Abstract 

This paper develops the small strain continuum dislocation theory accounting 
for statistically stored dislocations and Taylor hardening for single crystals. 
As illustration, the problem of anti-plane constrained shear of single crystal 
deforming in single slip is solved within the proposed theory. The distribution 
of geometrically necessary dislocations in the final state of equilibrium as 
well as the stress-strain curve exhibiting the Bauschinger translational work 
hardening and the size effect are found. Comparison with the stress-strain 
curve obtained from the continuum dislocation theory without statistically 
stored dislocations and Taylor hardening is provided. 

Keywords: dislocations, yield stress, hardening, stress-strain curve, size 
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1. Introduction 

Macroscopically observable plastic deformations in single crystals and 
poly crystalline materials are caused by nucleation, multiplication and motion 
of geometrically necessary dislocations (GNDs). There are various reasonable 
experimental evidences supporting the so-called low energy dislocation struc¬ 
ture (LEDS) hypothesis formulated first by Hansen and Kuhlmann-Wilsdorf 
rai- dislocations appear in the crystal lattice to reduce its energy (see also 
[261121]). Motion of dislocations yields the dissipation of energy which, in 
turn, results in a resistance to the dislocation motion. The general structure 
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of continuum dislocation theory (CDT) must therefore reflect this physical 
reality: energy decrease by nucleation of GNDs and resistance to the motion 
of GNDs due to dissipation. Just in recent years various phenomenological 
models of crystals with continuously distributed dislocations which are able 
to predict the density of GNDs as well as the accompanying size effects have 
been proposed in [T], [TTJ [T21 El El 0 EH E21 ESI HB1 UHl EOl [301 HS1 ESI ES] (see 
also the finite strain CDT proposed by Le and Stumpf [34,135], [36], Ortiz and 
Repetto [38j, Ortiz et al. [39], Le and Gunther [27], Koster et al. [21 ). 

In addition to the geometrically necessary dislocations there exists an¬ 
other family of dislocations which does not show up in the macroscopically 
observable plastic slip but nevertheless may have significant influences on 
the nucleation of GNDs and the work hardening of crystals. For any closed 
circuit surrounding an area, which is regarded as infinitesimal compared with 
the characteristic size of the macroscopic body but may still contains a large 
number of dislocations, the resultant Burgers vector of these dislocations 
always vanishes, so the closure failure caused by the incompatible plastic 
slip is not affected by them. Following Ashby [0] we call these dislocations 
statistically stored dislocations (SSDs). As a rule, the statistically stored 
dislocations in unloaded crystals at low temperatures exist in form of dis¬ 
location dipoles in two-dimensional case or small planar dislocation loops 
whose size is comparable with the atomic distance in three-dimensional case. 
The simple reason for this is that the energy of a dislocation dipole (or a small 
planar dislocation loop) is much smaller than that of dislocations apart, so 
the bounded state of dislocations renders low energy to the whole crystal. 
From the other side, due to their low energy, the dislocation dipoles (loops) 
can easily be created (as well as annihilated) by thermal fluctuations. The 
statistically stored dislocations play two important roles in the plastic defor¬ 
mations of crystals: i) together with the Frank-Read source p~5j they provide 
additional sources for the nucleation of GNDs due to the fact that, when the 
applied shear stress becomes large enough, the dislocation dipoles dissolve 
to form the freely moving GNDs, ii) the neutral dipoles (loops) of SSDs act 
as obstacles that impede the motion of GNDs leading to the nonlinear work 
hardening. In view of their important roles in ductile crystals, the account of 
SSDs in the CDT would make the material models more realistic. Arsenlis 
et al. [3] proposed a set of evolution equations for the densities of GNDs 
and SSDs within the crystal plasticity. The density of SSDs evolves through 
Burgers vector-conserving reactions, while that of GNDs evolves due to the 
divergence of dislocation fluxes. Except the missing thermal fluctuations in 
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the nucleation of SSDs, it was also unclear whether such an approach could 
be related to the energetics of crystals containing dislocations and the LEDS- 
hypothesis mentioned above. Berdichevsky [6j was the Erst who included the 
density of SSDs in the free energy density of the crystal. However, to the 
best of our knowledge, the more pronounced influence of the SSDs on the 
yield stress and the dissipation within the CDT has not been considered up 
to now. This paper aims at filling this gap. Its main idea is to propose the 
dissipation function depending on both densities of GNDs and SSDs in such 
a way that the obtained yield stress combines the constant plastic yield stress 
due to the Peierls barrier and the Taylor contribution that is proportional 
to the square root of the total density of GNDs and SSDs [12]. Then we 
apply the proposed theory to the problem of anti-plane constrained shear. 
We solve this problem numerically and find the distribution of GNDs in the 
final state of equilibrium as well as the stress-strain curve. We show the size 
effect for the threshold stress, the nonlinear work hardening due to the com¬ 
bined GNDs and SSDs, and the Bauschinger effect for the loading, elastic 
unloading, and loading in the opposite direction. 

The paper is organized as follows. In Section 2 the kinematics of CDT 
taking into account GNDs and SSDs is laid down. Section 3 proposes the 
thermodynamic framework for the CDT with SSDs and Taylor hardening. 
In Section 4 the problem of anti-plane constrained shear is analyzed. Section 
5 presents the numerical solution of this problem and discusses the distri¬ 
bution of GNDs, the stress-strain curve, the Bauschinger translational work 
hardening and the size effect. Finally, Section 6 concludes the paper. 

2. Kinematics 

In this paper we restrict ourselves to the small strain (or geometrically 
linear) continuum dislocation theory for single crystals. For simplicity we 
shall use some fixed rectangular cartesian coordinates and denote by x the 
position vector of a generic material point of the crystal. Kinematic quan¬ 
tities characterizing the observable deformation of this single crystal are the 
displacement held u(x) and the plastic distortion held /3(x) that is incom¬ 
patible. For single crystals having n active slip systems, the plastic distortion 
is in general given by 


|8(x) = V /3“(x)s“ ® m" (ft, = V (1) 

Cl=l d=l 
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with (3 a being the plastic slip, where the pair of constant and mutually or¬ 
thogonal unit vectors s° and m a is used to denote the slip direction and the 
normal to the slip planes of the corresponding a-th slip system, respectively. 
Thus, there are altogether 3 + n degrees of freedom at each point of this gen¬ 
eralized continuum. Here and later, equivalent formulas for the components 
of tensors are also given on the same line in brackets, where the Latin lower 
indices running from 1 to 3 indicate the projections onto the corresponding 
coordinates, while the Gothic upper index a running from 1 to n numer¬ 
ates the slip systems. We use Einstein’s summation convention, according 
to which summation over a repeated Latin index from 1 to 3 is understood. 
We see immediately from ([l]) that tr/3 = fin — 0, so the plastic distortion is 
volume preserving. 

The total compatible strain tensor field can be obtained from the dis¬ 
placement held according to 

£ = I(uV + Vu) (eij = ~ (Uij + U j:i )). (2) 

The incompatible plastic strain tensor field is the symmetric part of the 
plastic distortion field 

^=i(/3 + /3 T ) + &<))• P) 


For small strains we shall use the additive decomposition of the total strain 
into the elastic and plastic parts. Therefore, the elastic strain tensor field is 
equal to 


s e = 




(4) 


Although not absolutely necessary, we introduce for the illustration purpose 
the elastic distortion tensor field according to 


/3 e = uV-/3 (4;, //,, - d, ; ). 

The relationship between these three distortion fields is illustrated in Fig. [lj 
where F = I + uV, F p = I + /3, F c = I + (3 e . Looking at this Figure we 
see that the plastic distortion (3 is the distortion creating dislocations (either 
inside or at the boundary of the volume element) or changing their positions 
in the crystal without deforming the crystal lattice. In contrary, the elastic 
distortion /3 e deforms the crystal lattice having frozen dislocations [27]. 
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Figure 1: Additive decomposition 


Since the lattice rotation can directly be measured by the Electron Back- 
scatter Diffraction (EBSD) technique [25], it is interesting to express it in 
terms of the displacement gradient and plastic slip. Using the additive de¬ 
composition for the displacement gradients uV = e + uy we find the compo¬ 
nents of the total rotation tensor in the form 

w ij = 2 ( u hj — 

This tensor is obviously skew-symmetric. The plastic rotation tensor is the 
skew-symmetric part of the plastic distortion 

1 n 1 

<*$ = yA - ft*) = E 2^ s >i ~ s >*)' 

a=l 

The elastic (lattice) rotation tensor is the difference between these two tensors 


“ij = Uij 


uij = ^ u i,j 


u 




1 

V -p - »>?). 

Cl=l 


To characterize the geometrically necessary dislocations (GNDs) belong¬ 
ing to one slip system let us consider one term (3 a = /3 a s a <g> m n in the sum 
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0- The measure of incompatibility of this plastic distortion, introduced by 
Nye [37], Bilby [TO], and Kroner [22], reads 

« a = -/3“ x V (a?- = e jkl f3^ k ), 

where e l]k is the permutation symbol. With (3 a = (3 a s a © m“ we have 

«° = s° © (V/T x m a ) (a“- = s$e jk iP%mt). (5) 

To hnd the number of dislocations per unit area of the a-th slip system we 
take an infinitesimal area da with the unit normal vector n. Then the resul¬ 
tant Burgers vector of those geometrically necessary dislocations belonging 
to this slip system, whose dislocation lines cross the area da is given by 

b a = a" ■ nda (6“ = a^rij da). 

This is quite similar to the Cauchy formula relating the traction with the 
stress tensor. With cC = s a © (V/3 a x m a ) we get 

b a = s a ((V/3 a x m“) ■ n) da = s n (V/3 a • (m a x n)) da. 

Provided the direction t a tangential to the dislocation lines lying in the slip 
planes is known for the a-th slip system, then one can choose the infinitesimal 
area da with the unit normal vector t“ to compute the resultant Burgers 
vector of those geometrically necessary dislocations, whose dislocation lines 
cross this area under the right angle 

b“ = s“(V/P • (m a x t 0 )) da = s a d u (3 a da , 

where d u /3 a = V/3 a ■ v a , with is a = m n x t a being the vector lying in the slip 
plane and perpendicular to t“. This resultant Burgers vector can be decom¬ 
posed into the sum of two Burgers vectors: one parallel to t a representing 
the screw dislocations, and another perpendicular to t° corresponding to the 
edge dislocations. Therefore, the number of dislocations per unit area of each 
sort is given by 

= \ |s B - tfl ||0 I/ £°|, Py = ^l s °' v a \\dvP a \, (6) 

where the symbol © indicates the screw dislocations, while || the edge dislo¬ 
cations. 
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In addition to the geometrically necessary dislocations there exist another 
family of dislocations which does not show up in the macroscopically observ¬ 
able plastic distortion but nevertheless may have significant influences on 
the nucleation of GNDs and the work hardening of crystals. For any closed 
circuit surrounding an infinitesimal area (in the sense of continuum mechan¬ 
ics) the resultant Burgers vector of these dislocations always vanishes, so the 
closure failure caused by the incompatible plastic distortion is not affected 
by them. Following Ashby [Tj we call these dislocations statistically stored 
dislocations (SSDs). As a rule, the statistically stored dislocations in un¬ 
loaded crystals at low temperatures exist in form of dislocation dipoles (in 
two-dimensional case) or small planar dislocation loops, whose size is com¬ 
parable with the atomic distance (in three-dimensional case) (Arsenlis and 
Parks [2] have found also other three-dimensional self-terminating disloca¬ 
tion structures with zero net Burgers vector). The simple reason for this is 
that the energy of a dislocation dipole (or a small dislocation loops) is much 
smaller than that of dislocations apart, so this bounded state of dislocations 
renders low energy to the whole crystal. From the other side, due to their low 
energy, the dislocation dipoles can easily be created (as well as annihilated) 
by thermal fluctuations. As can be shown by the statistical mechanics of 
dislocations in two-dimensional case [8], the number of such dipoles remains 
nearly constant at the constant temperature. Let us denote the density of 
SSDs for each slip system by p a st . 

3. Thermodynamic framework 

To set up phenomenological models of crystals with continuously dis¬ 
tributed dislocation using the methods of non-equilibrium thermodynamics 
for irreversible processes let us begin with the free energy density. As a 
function of the state, the free energy density may depend only on the state 
variables. Following Kroner [23] we will assume that the elastic strain s e , 
the densities of GNDs a" (a = 1,... ,n), and the absolute temperature T 
characterize the current state of the crystal, so these quantities are the state 
variables of the continuum dislocation theory. The reason why the plastic 
distortion /3 cannot be qualified for the state variable is that it depends on 
the cut surfaces and consequently on the whole history of creating disloca¬ 
tions (for instance, climb or glide dislocations are created quite differently). 
Likewise, the gradient of plastic strain tensor e p cannot be used as the state 
variable by the same reason. In contrary, the dislocation densities ot a depend 
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only on the characteristics of GNDs in the current state (Burgers vector and 
positions of dislocation lines) and not on how they are created, so ad is the 
proper state variable. In addition to these state variable one should include 
also the densities of statistically stored dislocations p a st into the list of state 
variables. However, at low temperature these statistically stored dislocations 
prefer to exist in form of dislocation dipoles to render the crystal a low en¬ 
ergy. Provided the density of such dipoles depends only on the temperature, 
their energy contribution is a constant that can be omitted. Besides, due to 
the charge neutrality of dipoles, the energy of interaction between SSDs and 
GNDs are negligibly small compared to the energy contributions of GNDs. 
Thus, if we consider isothermal processes of deformation, then the free en¬ 
ergy per unit volume of crystal (assumed as macroscopically homogeneous) 
must be a function of s e and a“ 

ip = 'if{e e ,cx a ). (7) 

With this free energy density we can now write down the energy functional 
of the crystal. Let the undeformed single crystal occupy the region V of the 
three-dimensional euclidean point space. The boundary of this region, <9V, 
is assumed to be the closure of union of two non-intersecting surfaces, dk 
and d s . Let the displacement vector u(x) be a given smooth function of 
coordinates, and, consequently, the plastic slips /3“(x) vanish 

u(x) = u(x), /3 a (x) =0 for x G dk . (8) 

At the remaining part d s of the boundary the surface load (traction) f is 
specified. If no body force acts on this crystal, then its energy functional is 
defined as 

7[u(x),/5°(x)] = f ip(£ e ,ct a ) dx — / f u da, (9) 

JV Jd 3 

with dx = dxidx 2 dxz denoting the volume element and da the area element. 
Provided the resistance to the dislocation motion can be neglected, then 
the following variational principle turns out to be valid for single crystals: 
the true displacement held u(x) and the true plastic slips /3°(x) in the final 
state of deformation in equilibrium minimize energy functional ([9]) among all 
admissible fields satisfying the constraints (|8]). 

Let us find out the necessary conditions which must be satisfied by the 
energy minimizer in the final state of deformation in equilibrium. Taking the 


first variation of functional (J9]) we have 

51 = f [(TijiSeij - 5e^) + ^2 K ij^ a ij] dx - f fiSui da, 


where 


di\) 


° ij del 


Cl=l 


a = d±_ 

ij dot% 


a = 1 ,n. 


( 10 ) 


-y IJ 

We call cr (Cauchy) stress tensor, while n a higher order stress tensors. Taking 
the symmetry of a into account and using the kinematic relations ([l])-(|5|, 
we obtain 

P n n p 

51 = (cTijSuij - V + V' s*e jk i8fd\mf) dx - fi5u { da. 

J v n=1 a =i ' ^ 

Integration by parts with the use of the kinematic boundary conditions 
causing 5ui = 0 and 5f3 a = 0 on dk yields 


51= + s a i K a ij)k ei jk m a l )5[5 a ]dx 

J v o=l 

P n 

+ / [i a ij n j ~ fi)$ u i + 21. s i K ij e ijknkmi5/3 a ] da, 


a=l 


where r a = s^aijin^ is the resolved shear stress (Schmid stress). Since the 
variations 5ui and 5fd a can be chosen arbitrarily in V and on d s , from 51 = 0 
follow the equilibrium equations 

a ijd = 0, si^ewml + r a = 0, a = l,...,n (11) 

and the boundary conditions at d s 

o-ijUj = fi, s“K.“e /jfc n fc m / a = 0, a = 1 (12) 

The first equation of ( |TT| corresponds to the equilibrium of macro-forces act¬ 
ing on the volume element of crystal, while n remaining equations express 
the equilibrium of micro-forces acting on dislocations of the corresponding 
slip system. The first terms in ([II ]) 2 are called the back-stresses which are 
nothing else but the resultant forces acting on a dislocation from other dis¬ 
locations. Substituting (10) into ( [IT] ), we get the closed system of 3 + n 
governing equations to determine 3 + n unknown functions. 
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In real crystals there is however always the resistance to the dislocation 
motion causing the energy dissipation that changes the above variational 
principle as well as the equilibrium conditions. Various factors like impurities, 
inclusions, grain boundaries et cetera may have influence on this resistance 
to dislocation motion in polycrystals. However, in single crystals we may 
count two main sources. The first one is the periodic energy landscape of 
discrete crystal lattice (Peierls barriers) that, in combination with the small 
viscosity, leads to the rate independent dissipation 05- The second one is 
dislocations themselves acting as obstacles (dislocation forest) leading to the 
Taylor hardening |I2] . So the dissipation potential can be proposed in the 
form 

D = D(p%,pl,p°, t ,h- (13) 

We assume that this dissipation potential is the homogeneous function of the 
first order with respect to the plastic slip rates $ a . The simplest version of 
this dissipation potential would be 

n 

D = Y,S°(P^PpP’. t )\n (14) 

Cl=l 

where g a (p%, pjj, p%) are positive functions of the dislocation densities. In 
this case the latent hardening due to the cross-slip is neglected. The more 
complicated model taking the latent hardening into account is 


a,b=l 


where g ab (pQ, pfj, p a st ) is the n x n positive definite hardening matrix. 

When the dissipation is taken into account, the above formulated varia¬ 
tional principle must be modified. Following Sedov [33], Berdichevsky and 
Sedov [5] we require that the true displacement held u(x) and the true plastic 
slips /5“(x) in the final state of deformation in equilibrium obey the varia¬ 
tional equation 


51 + 



(15) 


for all variations of admissible fields u(x) and /3°(x) satisfying the constraints 
(|8]). Together with the above formula for 51 and the arbitrariness of 5ui and 
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8j3 a in V and on d s , equation (15) yields 


o. 


v,j 


= 0 , s a i K* jk e ljk m < l + T a = 


dD 

w* 


a = 1, 


,n 


(16) 


which are subjected to the boundary conditions ([ 8 ]) and (12). For the dissi¬ 
pation potential from (14) equations (16)2 become 

s i K ij,k € mmi +r a = g a {p a Q , pjj, /4) sign/) 0 , a = 1,..., n. (17) 

According to ([l7| the plastic slip /3 a may evolve only if the yield condition 


+ r “l = 0°(P©> Pp Pit) ( 18 ) 

is fulfilled. If |s°K 0 -^e;j fe m° + t°| < g a (pQ, pjj, p“ t ) then (3 a is “frozen”: /3 a = 0. 

In the following we shall consider single crystals deforming in single slip 
and take the simplest expressions for the free energy density mm 

-0 = ^A(tre 6 ) 2 + /itr(e e ■ s e ) + pk In -———, (19) 

2 1 - p/p s 


and dissipation 

D = (.K + g,ab\J p + p st )\$ \. (20) 

In these formulas p and A are Lame elastic constants, p s the saturated dis¬ 
location density, k and a are positive constants. It is assumed that only one 
sort of dislocation appears such that either p e = p (or p\\ = p). The second 
summand on the right-hand side of (20) describes the Taylor hardening due 
to the GNDs and SSDs. 


4. Anti-plane constrained shear 

As an application of the proposed theory let us consider the single crystal 
layer undergoing an anti-plane shear deformation. Let C be the cross section 
of the layer by planes z = const. For simplicity, we take C as a rectangle of 
width a and height h, 0 < x < a, 0 < y < h. We place this single crystal in a 
“hard” device with the prescribed displacement at the boundary dC x [0, L] 
(see Fig. § 

w = 'yy at dC x [0, L], 

where w(x,y,z ) is the ^-component of the displacement and 7 corresponds 
to the overall shear strain. We may imagine this single crystal as a grain, 
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Figure 2: Anti-plane constrained shear 


where the hard device models the grain boundary. The height of the cross 
section, h, and the length of the beam, L , are assumed to be much larger 
than the width a (a C ft, a C t) to neglect the end effects and to have the 
stresses and strains depending only on one variable x in the central part of 
the beam. If the shear strain is sufficiently small and the crystal is initially 
free of GNDs, then the crystal deforms elastically and w = 77/ everywhere 
in the specimen. If 7 exceeds some critical value, then the screw GNDs may 
appear. We admit only one active slip system, with the slip planes parallel 
to the plane y = 0 and the dislocation lines parallel to the z-axis. We aim at 
determining the distribution of GNDs as function of 7 within the framework 
of continuum dislocation theory taking into account the SSDs and Taylor 
hardening. For screw dislocations with the slip planes parallel to the plane 
y = 0, the tensor of plastic distortion, fiij, has only one non-zero component 
j3 zy = (3. We assume that (3 depends only on x-coordiriate: f3 = (3(x). Since 
the displacements are prescribed at the boundary of the crystal, dislocations 
cannot penetrate the boundaries x = 0 and x = a, therefore 

P(0) = /3(a) = 0. (21) 

The plastic strains are given by 

4 ? = 4 ? = 
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The only non-zero component of the tensor of dislocation density, o tJ = 

£jklPil,ki is 

&zz P,xi 


so the number of GNDs per unit area perpendicular to the z-axis is p = 
\P iX \/b. Under the assumptions made, the energy density per unit volume of 
the crystal with dislocations (19) takes a simple form 


U = -/i( 7 -/ 3) 2 + /ifcln 


1 - | A*I /Psb' 


( 22 ) 


while the dissipation potential becomes 


D 


IC + imbJ ,|/6 + p, t \f>\. 


(23) 


We first analyze the dislocation nucleation and dislocation pile-up during 
the loading when 7 increases from zero so that (3 > 0. In this case P{x) must 
satisfy the variational equation 


dx 


- PY + pk In -- 1 

to L 2 1 - \P,x\/Psb 

+ nabiJ\f3 tX \/b + p s ^j 5/3 dx = 0. 


(24) 


It is convenient to divide equation (24) by /i and introduce the following 
dimensionless quantities 


x = xbp s , a = abp s , 70 = —, £ = ctby/ps, k = —. (25) 

P Ps 

Since we shall deal with these dimensionless quantities only, the bar over 
them can be dropped for short. Equation (24) reduces to 

1 


2 ( 7 — PY + fcln — 


IAI 


dx+ ^70 + £ VW\ + K ) &P dx = 0 , 


( 28 ) 

with prime denoting the derivative. 

Due to the boundary conditions ( 21 ) P'(x) should change its sign on the 
interval ( 0 , a). One-dimensional theory of dislocation pile-ups [9] HSl J19,13TJ, 
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1521 35] suggests to seek the solution in the form 


A 0*0 

for x G (0, l), 



for x G (l, a — l), 

(27) 

Pi (a - x) 

for x G (a — l, a), 



where (3 m is a constant, l an unknown length, 0 < l < a/2, and P\{1) = /3 m 
at x = l. We have to find fti(x) and the constants, (3 m and l. Since (3[ > 0 
for x G (0,/), the first integral in (26) with this Ansatz reads 


'0 L 






dx + 2 (7 -/^m) 2 (a - 2/). 


Likewise the last integral in (26) becomes 


7o 


+ £V Pi + dx + (7o + ^a/k)(o - 2/) 5/? m . 


Varying the hrst integral with respect to /3i(x) and substituting the result 
in from it the equation for 

7 - 7o - WPi + K ~ A + 


into (26) we obtain from it the equation for /3i(x) on the interval x G (0,/) 

kfi'{ 


(1-W 


= 0. 


(28) 


The variation with respect to /3 m and l yields the two additional boundary 
conditions at x — l 


with 


W)= 0 , 2k — ( 7 r — 0 m )(a — 21) — 0, 


7r = 7 - 7o - S\ V 


(29) 

(30) 


Condition (29 )i means that the dislocation density must be continuous. 

At the onset of dislocation nucleation the dimensionless density of GNDs 
/3 1 must be small compared to 1, so equation (28) can be replaced by 

lr ~ Pi + Wl = 0. (31) 


Together with the boundary conditions (21 )i and (29)i this yields 

x lx 

Pi(' x ) = 7r(l — cosh —= + tanh —= sinh —=), 0 < x < l. 

Vk Vk Vk ~ ~ 
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Consequently, (29 ) 2 gives the following transcendental equation to determine 
l in terms of the constants k and a 


k l 

f(l) = 21 + 2 — cosh —= = a. 

lr Vk 


(32) 


According to (27) / must lie in the segment [ 0 , a/ 2 ]. Since cosh (l/Vk) > 1, 
2 1 < a — 2k/^ r . Thus, equation (32) has no positive root if 7 r < 2 k/a. 
Returning to the original variables according to (25) and (30) we see that 


inequality < 2 k/a corresponds to the condition 7 < 7 C , where 

2k 

lc = 7o + £\Ah-, 

a 

and for 7 < 7 C no dislocations are nucleated. This formula resembles the 
combined Taylor [42J and Hall [Pj], Petch {4U] relations. 

Assume now that at some stage of loading after the nucleation and ac¬ 
cumulation of GNDs we stop increasing 7 and then decrease it. The plastic 


slip rate could be either zero or negative. If the yield condition (18) cannot 


be fulfilled, then the plastic slip rate must be zero and the plastic slip f3(x) 
is frozen during this process. We call such process during which 7 decreases 
but j3(x) remains unchanged and equal to that plastic slip j3(x) at the end 
of the loading process elastic unloading. We will see that the average stress 
depends linearly on 7 during this clastic unloading. It remains now to ana¬ 
lyze the last case /3 < 0 which we call loading in the opposite direction. In 


this case the variational equation (26) changes to 


1 1 
W 


P) 2 + k In 


m 


dx 


(70 + £ VW\ + K ) ^ dx = 0 , (33) 


In the first stage of loading in the opposite direction the distribution of GNDs 
should not change much as compared to that at the end of the loading process. 
Therefore we can still assume /3(x) in the form (27). With this solution 


Ansatz it is straightforward to derive from (33) the following equation 


7 + 7o + £ VPi + K ~~ Pi + 


kP’{ 


(1 -P'i) 2 


= 0 . 


(34) 


and the boundary conditions 


Pi(l) = 0, 2k — (7r~ Pm)(a = 21), 


(35) 


15 















where 


7 1 = 7 + 7o + £a/«- 

If the dimensionless density of GNDs (3' is small compared to 1 at the end 


of the loading in the opposite direction, equation (34) can approximately be 
replaced by 

7 i~Pi + W{ = 0. (36) 


In this linearized version the system of equation (36) and boundary condi¬ 


tions (35) are identical with the system (31) and (29) if 7 1 is replaced by 
y r . Consequently, the solutions of these problems are equal if 7 1 = 7 r . In 
particular, the GNDs are completely annihilated at 

, r 2 k 

7 = -7c = -7o - $VK-• 

a 


5. Numerical simulations and comparison 



Figure 3: Evolution of /?: a) 7 = 0.005, b) 7 = 0.02, c) 7 = 0.05, d) 7 = 0.05 


In general the two-point boundary-value problem (28), (21)!, and (29)i 
for any given positive l G (0, a/2), due to its nonlinearity, can only be solved 
numerically. Since the slope of /3(x) at x = 0 is a priori unknown, we use the 
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Figure 4: Evolution of /3[(x) in the left boundary layer: a) 7 = 0.005, b) 7 = 0.02, c) 
7 = 0.05, d) 7 = 0.05 


collocation method and employ the Matlab program bvp4c for the solution 
of this two-point boundary-value problem [ni m- The third term of the 
equation (28) could lead sometimes to the complex solution if j3[+K suggested 
by the collocation method during the iteration turns out to be negative. 
Therefore we modify it to be + K ■ For 7 close to the threshold 

value where (3'(x ) is small, we use the solution of the linearized problem 
found above as the initial guess for the solution of the nonlinear equation. 
For larger 7 we solve the problem incrementally and use the solution obtained 
at the previous step as the initial guess for the solution at the next 7 differing 
from the previous one by a small increment. With this numerical solution 
at hand we can determine the function /3 m (l ) = so that equation (29h 
can be solved numerically. The latter is solved by the bisection method 
which is robust because it does not require the derivative of function on the 
left-hand side of (29 ) 2 - Fig. [3] shows the evolution of the plastic slip (3(x ) 
at different values of the overall shear 7. For the numerical simulation we 
took k = 10 -4 , p s = 10 15 m -2 , b = 10 _10 m, a = 10 _4 m, so that a = abp s = 
10. Other material parameters are chosen such that £ = 10~ 3 , k = 10 -2 , 
and 70 = K/p = 10 -3 . From this Figure one can see that the plastic slip 
/ 3(x ) increases as 7 increases. There are two boundary layers in which the 
geometrically necessary dislocations pile up against the obstacles at x = 0 
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and x = a. In the middle (/, a — l) the crystal is free of GNDs. 

a shows the evolution of the dimensionless density of geometrically 
necessary dislocations (GNDs) / 3 [(x) in the left boundary layer as 7 increases. 
We see that the density of GNDs is a monotonously decreasing function with 
the maximum being achieved at x = 0 . As 7 increases, the number of GNDs 
also increase and they pile up against the obstacle at x — 0. At 7 = 0.1 
the density of GNDs is already close to the saturated dislocation density. 
The dependence of the thickness of boundary layer l (where the GNDs are 
distributed) on the overall shear 7 is shown in Fig. | 5 } The thickness l starts 
from zero at 7 = 7 C and increases monotonically as 7 increases. For moderate 
and large values of 7 the thickness / is nearly a linear function of 7. 


I 



Figure 5: Thickness of the boundary layer l as function of 7 


As soon as the plastic slip develops, the shear stress a = //(7 — ( 3 (x)) 
becomes inhomogeneous. It is interesting to plot the average shear stress 

1 r a 

a = - /x(7 — / 3 (x)) dx ( 37 ) 

a Jo 


as function of the shear strain. For 7 < 7 c the plastic slip / 3 (x) = 0 , so 
a = fi'f. For 7 > 7 c the average stress will be less than /i 7 due to the 
positiveness of the plastic slip. However, one can still observe the positive 
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slope of the stress-strain curve which can be explained by the kinematic and 
Taylor hardening in combination. Fig. [6] shows the normalized average shear 
stress cr/p versus the shear strain curve (bold line). As comparison, we show 
also the stress-strain curve obtained in the continuum dislocation theory 
without SSDs and Taylor hardening (dashed line) as well as the linearized 
theory without SSDs and Taylor hardening (dashed and dotted line). One can 
see that the slope of the stress-strain curve obtained by the CDT taking into 
account SSDs and Taylor hardening is highest. Mention that the threshold 
stress at which the GNDs begin to nucleate in the theory without SSDs and 
Taylor hardening is 70 + 2k/a that is also lower than the threshold value 7 C . 

0/|X 



Figure 6: Normalized average stress versus shear strain curve: i) theory with SSDs and 
Taylor hardening (bold line), ii) theory without SSDs and Taylor hardening (dashed line), 
iii) linearized theory without SSDs and Taylor hardening (dashed and dotted line) 


On Fig. [7] we show the normalized stress-strain curve obtained during 
loading, elastic loading and unloading, and loading in the opposite direc¬ 
tion. The curve AB corresponds to the stress-strain curve during the loading 
process when 7 increases from 7 C . The straight lines DA going through the 
origin and BC are the stress-strain curves during elastic loading and unload¬ 
ing when 7 increases (decreases) at frozen /3. It can be seen from (37) that 
the slope of DA and BC must be constant because of the frozen /3(x). Finally, 
the curve CD corresponds to the stress-strain curve during the loading in the 
opposite direction when the plastic slip /3(x ) decreases simultaneously with 
the decrease of 7. This curve is obtained from (37) where the plastic slip 
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(3(x) corresponds to the solution of the problem (34) and (35). The latter is 
solved numerically with the use of Matlab bvp4c in the similar manner as for 
the system (28) and (29). This stress strain curve cuts the elastic line AD at 
point D with the coordinate (—y c , — 7 C ) at which the plastic slip j3(x) must 
be equal to zero. 



Figure 7: The stress-strain curve during loading, elastic loading and unloading, and loading 
in the opposite direction 

Let us consider now the following close loading path: 7 is first increased 
from zero to some value 7 * > 7 C , then decreased to — 7 C = —70 — — yb 

and finally increased to zero (Fig. [ 8 ]). The rate of change of 7 (t) does not 
affect the results due to the rate independence of the dissipation. 

y 



Figure 8: Loading path 

In Fig. [7] the close loop OABCDO shows the normalized average shear 
stress (or average elastic shear strain) versus shear strain curve for the load¬ 
ing program of Fig. [ 8 | The straight line OA corresponds to the purely elastic 
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loading with 7 increasing from zero to 7 C , where (3(x) = 0. The line AB 
corresponds to the plastic yielding where the yield condition (Il8|) (or, equiv¬ 


alently, the system (28) and (29)) holds true. The yield begins at the point 
A with the yield stress a y = p 7 c . The work hardening due to the nucleation 
and pile-up of GNDs and Taylor hardening is observed. During the unloading 
as 7 decreases from 7* to 7* = 7* — 2^ c (the line BC) the plastic slip (3 = (3* 
is frozen. As 7 decreases further from 7* to — j c , the plastic yielding occurs 
where the yield condition (18) (or, equivalently, the system (34) and (35)) 
holds true (the line CD). The yield stress a y = /i( 7 * — 2 7 J at the point C, at 
which the inverse plastic flow sets on, is larger than — /ry c (because 7* > %). 
Along the line CD, as 7 is decreased, the created GNDs annihilate, and at 
the point D all GNDs disappear. Finally, as 7 increases from — j c to zero, the 
crystal behaves elastically with [3 = 0. In this close cycle ABCD dissipation 
occurs only on the lines AB and CD. It is interesting that the lines DA and 
BC are parallel and have the same length. In phenomenological plasticity 
theory this property is modeled as the translational shift of the yield surface 
in the stress space, the so-called Bauschinger effect. 


6. Conclusion 

In this paper we have proposed the continuum dislocation theory taking 
into account the statistically stored dislocations and Taylor hardening. On 
the example of the anti-plane constrained shear we have shown that the 
threshold stress at which the geometrically necessary dislocations nucleate 
is higher than that of the theory without SSDs and Taylor hardening. For 
the loading path consisting of increasing and decreasing shear strain the 
stress-strain curve becomes a hysteresis loop. The work hardening and the 
Bauschinger effect are quantitatively described in terms of the plastic slip 
and the densities of GNDs and SSDs. 
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